Phase vortices from a Young's three-pinhole interferometer 
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An analysis is presented of the phase vortices generated in the far field, by an arbitrary arrangement of 
three monochromatic point sources of complex spherical waves. In contrast with the case of three interfering 
plane waves, in which an infinitely-extended vortex lattice is generated, the spherical sources generate a finite 
number of phase vortices. Analytical expressions for the vortex core locations are developed and shown to 
have a convenient representation in a discrete parameter space. Our analysis may be mapped onto the case 
of a coherently-illuminated Young's interferometer, in which the screen is punctured by three rather than two 
pinholes. 
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I. INTRODUCTION 

In a seminal paper, Dirac Q] considered vortical screw- 
type dislocations in the phase of complex wavefields, not- 
ing the one-dimensional nature of the associated vortex cores 
(nodal lines) in three dimensions (3D). Such phase vortices 
exist in a variety of linear and non-linear physical systems that 
may be described via complex fields, including the angular- 
momentum eigenstates of the hydrogen atom the Meiss- 
ner state of type-II superconductors iS 3], vortex states of 
superfluids (a, \S\ and Bose-Einstein condensates [7], opti- 
cal vortex solitons [8], propagating electron wavefunctions 
diffracting through crystalline slabs |2D, Gaussian random 
wavefields [10] and optical speckle fields Hill . 

In continuous complex scalar fields, to which the consid- 
erations of the present paper are restricted, vortical behav- 
ior is manifest as screw-type dislocations in the field's multi- 
valued surfaces of constant phase (13, E, [HI]. More pre- 
cisely, consider a stationary-state, complex spatial wavefunc- 
tion or order-parameter field ^(r) = A(r)exp [/^(r)]. Here, 
A(r) is the non-negative real amplitude, ^(r) is the phase, and 
r = (x, y, z) is a position vector in 3D. Note that harmonic time 
dependence on angular frequency a> and time t, via the usual 
multiplicative factor of exp(-ia)t), is suppressed throughout. 
To determine whether a phase vortex exists at a point p in 
a plane II over which ^(r) is defined, a line integral of the 
phase gradient is evaluated along a smooth infinitesimally- 
small path T that encircles p. This path is assumed to have 
a winding number of unity with respect to p, and to be such 
that the modulus of ^(r) is strictly positive at each point on T. 
One may then write the following expression for the "circula- 
tion" of the field over T (see, e.g., Nye lfl5ll ): 



i dx= i 



Vx - ids = 2m. 



(1) 



Here, t is a unit tangent vector to T, ds is an infinitesimal line 
element along T, and n is an integer corresponding to the num- 
ber of phase windings about p. Any non-zero n indicates the 
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presence of a vortex core threading the path T, with the non- 
zero value for n being referred to as its topological charge. 
The sign of this charge distinguishes between a vortex (+) and 
an anti- vortex (-). 

In an optical setting, a common way to generate such phase 
vortices is to pass coherent laser or soft X-rays through a 
spiral phase plat e or forked transmission diffraction grating 
E O EE M HI [H. Kim et al. HI generated opti- 
cal vortices with a curved glass plate as an alternative to the 
spiral phase plate. Another means for generating fields with 
specified vortical phase is a dynamically computer-controlled, 
holographic system ll23ll . Other means for creating optical 
phase vortices in coherent light include the use of spatial light 
modulators HI], the use of aberrated lenses to create vortices 
in a distorted focal volume i25h . and diffraction from random 
phase screens fToh . 

As an alternative approach, one can forego the use of 
diffractive or refractive optical elements, seeking instead to 
create phase vortices by the superposition of a small number 
of non- vortical fields. For example, Nicholls and Nye B26ll 
showed that one can generate a lattice of vortices by inter- 
fering three plane waves. Later, Masajada and Dubik i27ll 
showed that this is a minimum requirement and reformulated 
the analysis in terms of phasors. 

Here we generalize this idea, by considering the forma- 
tion of phase vortices via the superposition of three outgoing 
spherical waves, generated by three distinct monochromatic 
equal-energy point sources. We see that the resulting system 
of nodal lines (vortex cores) exhibits a rich geometry, by de- 
veloping approximate analytical expressions for the far-field 
behaviour of this nodal-line network. The three-dimensional 
space, into which the sources radiate, is foliated using a family 
of observation planes that are parallel to the plane containing 
the three point sources. When one observes the wavefield over 
any such foliating plane, a 2D pattern of point vortices may be 
seen, the cores of which coincide with the points at which a 
nodal line punctures the plane. 

Interestingly, the problem of three interfering spherical 
waves may be mapped onto a Young-type experiment, in 
which a black screen with three small pinholes is coherently 
illuminated by a propagating complex scalar field. Note that 
this identification is only possible when one is both suffi- 
ciently far from the screen and sufficiently close to the optic 
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axis, in which case the radiation transmitted by each of the 
pinholes is approximately spherical. 

We close this introduction with a brief outline of the re- 
mainder of the paper: We begin by reviewing the manner in 
which the superposition of three plane waves may be used 
to generate an infinite lattice of phase vortices. The gener- 
alization of this idea, to the superposition of three outgoing 
spherical waves, is then given. We describe the application 
of a phasor approach to the spherical- wave arrangement, ap- 
plying this in the far-field region asymptotically far from the 
sources. Approximate analytical expressions are derived for 
the vortex locations. A representation in terms of a certain 
parameter space arises, allowing estimates of the number of 
vortices and description of a natural coordinate system for the 
vortices at the intersections of a certain family of hyperbolas. 
A specific case of collinear sources is explored in detail. We 
then show how the theory, which has been derived for spheri- 
cal point sources, may be mapped onto the case of a Young's 
interferometer in which the illuminated screen contains three 
rather than two pinholes. 
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II. PHASE VORTICES FROM THE INTERFERENCE OF 
THREE PLANE WAVES 

Consider the following superposition of three planar spatial 
wavefunctions: 



Y(r) = 2^exp[/(k r r + y -)] > 



(2) 



7=1 



where the non-negative real constants Aj denote the ampli- 
tude of the 7th wave, k 7 are wavevectors corresponding to the 
same de Broglie wavelength Aq = 2n/\kj\, and (pj are global 
phase factors. Notwithstanding the fact that the constituent 
plane waves do not have a vortical character, the above su- 
perposition may yield a regular lattice of phase vortices and 
anti- vortices 11261 12711. 

A numerical example of this phenomenon is given in Fig.[TJ 
corresponding to the parameters A\ - A^ - A3 = 1 and 
01 = 02 = 03 = 0, with all fields being evaluated over 
the plane z = 0. This example illustrates the three interfer- 
ing plane waves giving rise to an infinitely extended lattice of 
straight, parallel nodal lines (vortex cores). These nodal lines 
intersect the plane z = 0, to give the location of the point 
vortices that are visible as screw dislocations in the phase 
map of Fig. Old). The locations of these point vortices co- 
incide with both (i) the intersections of the zero contours in 
Fig- Hlb), and (ii) the points at which |^| = in the amplitude 
plot of Fig. [He). Indeed, continuity of the wavef unction im- 
plies that the probability density vanishes at each vortex core, 
since these are branch points at which the phase ceases to be 
differentiable. The topological charge n of each of these point 
defects is seen to be equal to ± 1 , as the phase increases by ±2n 
as one traverses a circuit that encloses a given vortex core. In 
this context, we note that higher-charge vortices are unstable 
with respect to perturbation [28], which is why they are not 
observed in the present setting. 



FIG. 1: (Color online) Interference of three plane waves giving 
rise to an infinite, regular vortex lattice [see Eq. ©]. Here the 
wavevectors are oriented symmetrically with respect to the z axis, 
with ki = (V3,-l,10),k 2 = (- V3, -1, 10) and k 3 = (0,2,10). 
The real (a) and imaginary (b) parts, amplitude (c), and phase (d) 
of the wavefunction ^(x, y,z = 0) are shown. In (b), zero contours of 
the real (dashed lines) and imaginary (solid lines) parts are overlaid. 
In (c) is overlaid the vector field of the probability current density, 
which is seen to rotate counter-clockwise around vortices and clock- 
wise around anti- vortices. In (d), vortices are visible as the ends of 
branch cuts indicated by dark to light (-71 — > n) steps. Values in 
(a-c) are represented by linear levels from dark to light (minimum to 
maximum). 



III. PHASE VORTICES FROM THE INTERFERENCE OF 
THREE SPHERICAL WAVES 

Given that the superposition of three complex plane-wave 
spatial wavefunctions may lead to phase vortices J26[ |27|], it 
is natural to enquire whether the superposition of three outgo- 
ing spherical waves may not also lead to phase vortices. This 
latter case is investigated here. 



A. Extending the plane- wave case to spherical waves 

The complex spatial wavefunction ^(r), due to a point 
source at position r 7 - that is radiating outgoing spherical waves 
in vacuo, is given by 



W=7^exp[/(^|r-r 7 | + 7 )] 



—. exp (ix) , 



(3) 



r-r,- 



where k = In/A®. For all r ± r 7 , such spherical waves 
obey a variety of linear partial differential equations, in- 



FIG. 2: Coordinate system. The point sources r\, r 2 , and r 3 lie in the 
plane z = 0. The position vector r is to a point with cylindrical-polar 
coordinates (r ± ,6, zo)- 



eluding: (i) the time-independent free-space Schrodinger 
equation for non-relativistic spinless particles, (ii) the time- 
independent free-space Klein-Gordon equation for relativistic 
spinless particles, and (iii) the free-space Helmholtz equation 
for monochromatic complex scalar electromagnetic waves. 
As such, the following discussions are applicable to all of 
these physical systems. 

Now consider an assembly of three point sources, all of 
which have the same wave-number k. Without loss of gener- 
ality we may consider these sources to occupy the same plane 
z = 0, with source locations r ; = (xj,yj, 0), where j - 1, 2, 3. 
The resulting spatial wavefunction ^(r) may thus be written 
as 



j_j | r _ r .| L V 



■r/ + 



(4) 



where |r — r 7 | = ^{x- Xj) 2 + ( y - yj) 2 + z 2 is the distance 
from the jth source to a given observation point r. 

Referring to Fig. O we label both the 7th source and its 
distance from the coordinate origin by the same symbol r 7 . 
The position vector r and its perpendicular component r ± have 
lengths r = |r| and r± = \r±\, respectively. 

To determine the location of the vortices which result from 
the superposition of the three spherical waves, we utilize the 
fact that vortex cores lie at points of zero amplitude. Given 
that the problem is restricted to two degrees of freedom, due to 
the complex wavefield representation, a geometric phasor dia- 
gram can be constructed with one phasor for each wave com- 
ponent in Eq. ©: see Fig. [3] We follow the phasor approach 
of Masajada and Dubik [27] (see also J29D). At any point co- 
inciding with a vortex core, the phasor components must sum 
to zero when placed tip to tail. Note that if the source ampli- 
tudes differ sufficiently, it is possible that no closed triangle of 
phasors may be formed, i.e., circles A and B in Fig.[3]cannot 
intersect if |pj > |p 2 | + |p 3 |. In this case, no vortices will be 
produced. 

Consider a given point r in space, corresponding to the spe- 
cial case of Eq. © where A\ = A2 = A3 = 1. Let p 1? p 2 ,p 3 
denote the three complex terms that are summed on the right 
side of this equation. These three numbers are represented as 
phasors in Fig. [3] Here, we set (pi = 0, which implies no loss 
of generality, since the invariance of the equations of motion 
under a shift in the origin of time implies global phase fac- 
tors to have no physical meaning. The circle B ' represents the 



FIG. 3: Phasor diagram from which vortex solution conditions are 
established [27]. Three phasors p 1? p 2 , p 3 — which respectively cor- 
respond to the three terms on the right side of Eq. ©, at a given point 
r on a nodal line — form a closed triangle corresponding to a resul- 
tant zero amplitude. Two angles y and 77 arise in the construction, as 
indicated, (a) In the case of arbitrary amplitudes for p 1? p 2 and p 3 , 
the tip of p 2 is constrained to lie on circle A of radius |p 2 |, with the 
tail of p 3 being constrained to lie on circle B of radius |p 3 |. (b) For 
equal amplitudes |p 1 1 = |p 2 | = |p 3 | an equilateral triangle is formed. 
The dashed construction represents another of the six equivalent al- 
ternatives formed by permutating the phasor order. 



possible orientations of p 3 , constrained by the tip of p 2 . The 
zero sum condition is easier to construct if p 3 is flipped and 
the circle B drawn constrained by the tail of \) l . The resulting 
vortex solutions correspond to the phasors p 2 and p 3 meeting 
at the intersections of the circles A and B. Note that there 
are two such "closed triangle" intersections and hence two 
apparent solutions. In the case of equal amplitudes, where 
IPil = IP2I = IP3I' symmetry dictates that there are only two 
unique solution angles. We consider the triangle in the first 
quadrant of the complex plane, which is an arbitrary choice, as 
any of the six equivalent constructions formed by permutating 
the phasor order will lead to the same solution. Relaxing the 
equal-amplitude condition will require consideration of extra 
solutions. 



B. Vortices in the far-field regime 

We evaluate Eq. © in the "far field" regime, namely in the 
half- space z > Zo in which zo is sufficiently large that 



r-r,- 



-, r = r » r 7 - 



(5) 



This approximation is applied to the phase term in the expo- 
nent of Eq. (|4]). The wavefunction *F depends linearly on the 
amplitude term Aj/\r - r 7 -|, so for large r the divisor varies 
much more slowly with |r - r 7 | than the phase argument. Con- 
sequently, the stronger approximation |r — r 7 | ~ r is made to 
this term (see, e.g., L30]). Thus, Eq. d4j) becomes 



Y(r) 



exp<z 



3 1 

Y- 

- exp (ikr) + - exp \ i 
r r 



+ - exp \ 



k\r cos 9 J + 

cos(#-0 3 )) + ^ 



(6) 



4 



where we have made use of the assumptions that A\ = A 2 = 
A 3 = 1 and that the sources share a single wavenumber k. This 
expression vanishes when 



1 + exp 
+ exp 



i (-k-^ cos + cp2 j 

/(-£^COS(0-0 3 )+0 3 ) 



(7) 



= 0. 



Geometrically, the above condition reduces to the addition 
of three unit-length phasors in the complex plane, such that 
they form an equilateral triangle when placed tip-to-tail l27ll . 
This construction is shown in Fig. [^b), with the arguments 
of the two exponentials in Eq. © being denoted by y and 
77, respectively. These phase angles are uniquely defined to 
within an integer multiple of 2n, so that: 



r ± r 2 

-k COS + 02 

r 



2tt 



+ 2mn 



and 



-k^- 1 cos (0 - 63) + 03 = — + 2nn, 
r 3 



(8a) 



(8b) 



where m and n are integers. 



C. Vortex locations 

Here, we show how the construction of the previous sub- 
section can be used to determine the polar coordinates r ±mn 
and mn , of a given point vortex in the plane z = zo, that is 
specified by the integer indices (m, n) (cf. Fig. [2]). 

Dividing Eq. (l8bl) by (l8at gives 

An 



r 3 cos(0-0 3 ) ^ y + 2 ^"^3 

r 2 cos0 27r 

— + 2mn - 02 

2(2 + 3n)n - 30 3 



(9) 



2(1 + 3m)7T-30 2 



We denote the denominator and numerator, on the right-hand 
side, as 



M(m) = 2(1 + 3m)/r- 302, 



and 



N(ri) = 2(2 + 3n)n - 30 3 
respectively. Next, making the substitution 

cos (0 - 63) cos cos 63 + sin sin 3 
cos cos 

= cos 3 + tan sin 3 , 



(10a) 



(10b) 



(11) 



gives 



r 3 I \ N(n) 

— cos0 3 +tan0sm0 3 = — — . 
r 2 \ / M{m) 



(12) 



Finally, isolating and labelling it with an mn subscript, to 
identify it with the (m, ft)th vortex core, gives the desired ex- 
pression for the polar angle to the (m, ft)th vortex core, 



6 mn = arctan 



1 



r 2 N(n) 
sin0 3 \r 3 Mini) 



■ cos 3 



(13) 



With a view to obtaining the radial coordinate r ±mn of the 
vortex core, take Eq. (l8at and write the denominator r in terms 
of its components zo and r ± (see Fig. [2]). Hence: 



r±r 2 



: COS m 



1 

-Mim). 



V^o 2 + r ± 2 

Squaring, and then solving for r ± , we obtain 

2 



^ _Lmn — 



^0 



3&r 2 cos 6 m 
M(m) 



(14) 



(15) 



- 1 



where an mn subscript has been added to r ± . Applying the 
identity 



cos mn — 



1 



1 + tan 2 6 m 



(16) 



and making use of Eq. (fT3l) , we obtain our final expression for 
the radial coordinate r ±mn of the (m, ^)th vortex core: 



(3kr 2 /M(m)Y 



(17) 



1 + 



1 



r 2 N(ri) 



sin 2 3 \r 3 M(m) 



- cos 3 



The positive and negative solutions correspond to two separate 
vortices. Note that these coincide with an extra branch of the 
arctan function in Eq. (TT3l) . Note, also, that Eq. (TTTb is only 
valid for integers (ra, n) that yield a real number for r ±mn (cf . 
Sec.HIlD]). 

The polar equations (IT3l) and (IT71) specify the vortex core 
locations for all allowed m and n parameter values, in the far- 
field regime. Note that r ±mw is proportional to zo, as one would 
expect in the far-field. This may be contrasted with the case 
of three superposedplane waves, where the nodal lines are 
mutually parallel I26i l27ll. 



D. Parameter Space 

For real solutions, the argument of the square root in 
Eq. dT71) must be positive, imposing a condition on the allow- 
able (m, n) values for a given source arrangement. In what fol- 
lows we set cpj = 0, corresponding to all three point sources 
radiating in phase with one another. The integers m and n must 
therefore satisfy the inequality 



( 1 + 3m) sin fh | + | (2 + 3n) y - (1 + 3m) cos 3 j 



3&r2 sin 3 N 
2n , 



(18) 



We claim that this describes the interior of an ellipse in the 
Cartesian (m, n) plane, for all non-collinear arrangements of 
the three sources. 

To prove the above claim, first note that the boundary curve, 
of the open region defined by Eq. (fT8l) . is obtained by replac- 
ing the inequality in this expression with an equality. The re- 
sulting equation is consistent with the form of a general conic 
section in the (m, n) plane, namely (see, e.g., 113111 ): 



Q(m, n) - am 2 + 2hmn + bn 2 + 2gm + 2fn + c = 0, 
where a,h,b,g,f and c are real numbers given by 
a = 9r 3 2 

h = -9r2r3 cos 63 
b = 9r 2 2 

g = 3r3 2 - 6r2r3 cos 63 
f = 6r2 2 - 3r2r3 cos 63 

9,9, /3&r2r3sin#3 
c = r 3 2 + 4r 2 2 - 4r 2 r 3 cos6> 3 - -j^ 



a h g 
h b f 
g f c 



(19) 



(20) 



Introduce the invariants Oil 



A = 



a h 
h b 



(21a) 



(21b) 



and 



r = a + b. 

Substituting Eqs (l20t into Eqs ([211 and evaluating gives 

^27& ,^ 2 
A = - I — (r 2 r 3 sin# 3 ) : 



(21c) 



and 



6 = (9r 2 r 3 sin6/3) 2 



r = 9(r 2 2 + r 3 2 ). 



(22c) 



In order for Eq. (IT9l) to correspond to an ellipse, the discrim- 
inant conditions A ^ 0, 5 > and A/r < must be satisfied. 
These three conditions are met when: (i) r 2 , r3 ± 0, and (ii) 
63 ± pn, p e Z. This will always be true for non-collinear 
arrangements of three distinct sources. Since the area of the 
corresponding ellipse is finite, for three non-collinear sources 
each of which are separated by a finite distance, we have a 
finite number of vortices labelled by the integer pairs (m, n) 
obeying Eq. (fT8l) . 

The parameter- space ellipse has center (mo, no), rotated 
anti-clockwise at an angle (p, with semi-axis lengths a' and 
b' (Fig.©. These ellipse parameters are given by well known 




FIG. 4: A parameter- space ellipse, in the (m, n) plane, sets an upper 
limit on the number of vortices that are created by the interfering ra- 
diation from three point sources, (a) The bounding ellipse has center 
(mo, wo), makes an angle <p to the positive m-axis and has semi-axis 
lengths a' and b' . (b) The enclosed (m,n) pairs lie on a discrete 
square lattice of unit spacing. Each of these interior points corre- 
sponds to two vortices. 



expressions in terms of the coefficients in Eq. (l20l) . The center 
and rotation angle are given by 113 lh 



(mo, no) = 



bg-hf af-hg 
h 2 - ab' h 2 - ab 



1 lb-a\ 1 

= 2 arccot (^r) = 2 



1 2 
~3'~3 



r 3 2 - r 2 2 

^ arccot I 

2 \2r2r3COS#3 



(23) 



The semi-axis lengths a' = s+ and b' = S- are given by 



A 



(24) 



(25) 



where A+ denotes the two solutions to the quadratic 
A 2 - tX + 6 = 0. 

Thus 

A± = \ (ri 2 + r 3 2 ± A /(2r2r 3 cos^) 2 + (r2 2 -r 3 2 ) 2 J , (26) 
(22a) so mat me sem i -ax i s lengths are given by 



s ± = 



|fc^2 r 3 sin 03 1 



(22b) n ^2 A /(2r 2 r3COS03) 2 + (r2 2 -r 3 2 ) 2 ± (r 2 2 + r 3 2 ) 



(27) 

For fixed k, r 2 , and ¥3 , and variable 63 , Eqs (l23t and (l27l) de- 
fine a family of ellipses. A "bounding rectangle" may be con- 
structed, as the envelope of this continuum of ellipses. Any 
one ellipse in this family, corresponding to a particular value 
of 63, touches each side of this bounding rectangle exactly 
once. The bounding rectangle is centered at (mo, no), and has 
dimensions of kr 2 /n and kr3/n in the m and n directions, re- 
spectively. (Note that these dimensions are found by setting 
63 = 71 12 in Eq. (127b.) As 63 is varied from to n, while 
keeping k, r 2 , and r3 fixed, the parameter- space ellipse trans- 
forms from: (i) a line at 45° to the m-axis, identified with the 
positive-gradient diagonal to the bounding rectangle, to (ii) a 
series of non-degenerate ellipses, each of which touch each 
side of the bounding rectangle exactly once, to (iii) a line at 
-45°, namely the negative-gradient diagonal to the bounding 




FIG. 5: (Color online) (a,c) Amplitude |*F| and (b,d) associated phase 
X from Eq. © at zq-25Aq with different source arrangements, shown 
pictographically alongside the (a) and (c) labels. Small circles over- 
laid on numerical simulations show the vortex locations as approxi- 
mated by Eqs (TTTb and fT71 ). Values in (a, c) are represented by linear 
levels from dark to light (minimum to maximum). The phase ranges 
from -71 (dark) to n (light). All sources are in phase (i.e. 02=03=0), 
Ao = 7i and the source arrangements are: (a,b) r 2 = r 3 = 3/t , # 3 = 0°. 
(c,d) r 2 =r3=3Ao, 03=60°. Note that a spherical background has been 
subtracted from all phase maps, as described in the main text. 



rectangle. Note that no vortices are produced in the limit cases 
(i) and (iii) above, since the open region bounded by a straight 
line is an empty set (cf. Eq. ([T8V). The change in shape of the 
ellipse with 63 is symmetric about 63=71 (corresponding to 
case (iii)), so that the ellipse for 63 = 71 - £ is coincident with 
that for 71 + for any angle 

Figs [5] and [6] present simulations with various source ge- 
ometries showing predicted vortex locations for the far-field 
case. The corresponding parameter space ellipses are shown 
in Fig. [71 The fields of view of the intensity and phase plots 
do not show the outermost vortex cores in most cases. This 
is evident from a count of the lattice points enclosed by the 
corresponding ellipse in Fig. [7] Only Figs 0c) and (d) have 
all the generated vortices within the visible region. The cor- 
responding parameter- space ellipse, near -45° in Fig. [71 en- 
closes 6 lattice points — this corresponds to the 6 x 2 = 12 
phase vortices in FigsOc) and (d). 

Note that, to aid visualization in all of the phase plots 
in Figs [5] and [6l a constant spherical background has been 
subtracted. Indeed, far from the three point sources, one 
may meaningfully write the wavefunction as a single expand- 
ing "background" spherical wave, multiplied by an envelope 
whose functional form depends on the particular local ar- 
rangement of the sources. By subtracting the phase of this 
spherical background from all of the displayed phase maps, 
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FIG. 6: (Color online) Phase x from Eq. $4$ and vortex-core loca- 
tions approximated by Eqs JT3l > and (Q7}, with variation in source 
arrangement. The representation and parameters are as described in 
Fig. [5] The scale-bar in (c) applies to all images and corresponds 
to zo = 25/lo, except for (d) which corresponds to zo = 250/io- Small 
arrows in (b,c) point to the exact vortex locations, to indicate their 
deviation from the far-field predictions. Aq=7t and source arrange- 
ments are: (a) r 2 = r 3 = 3A , 6 3 = 10°. (b) r 2 = r 3 = 3A , 63 = 90°. 
(c) r 2 = r 3 = 3A , 6» 3 = 175°. (d) r 2 = r 3 = 3A , 6> 3 = 175°, zo = 25(U . 
(e) r 2 =r 3 =3A , 6> 3 = 180°. (f) r 2 =3Ao, r 3 =1.5/lo, 6> 3 =60°. A spherical 
background has been subtracted from all phase maps (see main text). 



the structure of the envelope alone — including any vorti- 
cal structure — may be examined, without the distraction of 
a large number of concentric phase contours from the back- 
ground wave. Note also that this background spherical wave 
corresponds to an effective source located at the geometric 
centroid of the real sources, with amplitude A = A\ + A 2 + A3. 

Return consideration to Fig. [3a), which exhibits a degen- 
erate case in which two of the point sources are co-located, 
thereby reducing the system to two in-phase spherical sources 
with one having twice the amplitude of the other. As expected 
for only two sources, no vortices are generated. Rather, one 
has a series of Young-type fringes, with minima correspond- 
ing to a series of nodal planes in 3D. Each of these nodal 
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FIG. 7: Parameter space ellipses corresponding to Figs [5] and [6] En- 
closed lattice points define vortex locations. 



planes is seen to constitute a "domain wall" for the phase 
|0] [see Fig. Ob)]. The phase is not left-right mirror sym- 
metric, due to the positioning of the geometric centroid for 
three sources being to the right of center. Thus, there is a 
global tilt to the phase, which is observed to cycle through two 
branch-cuts. For this example, the parameter- space ellipse 
(not shown) is the limiting case of a line at 45°. The comple- 
mentary limiting case is Fig.Oe) where the three sources are 
collinear; the ellipse (not shown) is a line at -45° and again 
no vortices are generated. The phase here is mirror symmetric 
since: (i) the initial field configuration is mirror symmetric, 
and (ii) the wavefield propagator is rotationally symmetric. 
In Fig. Oa) one of the sources has been moved just enough 
from coincidence with another so that the corresponding el- 
lipse just encompasses some lattice points and vortices are 
created. There are four points clearly within the ellipse, giv- 
ing rise to the vortices seen in the panel. If the additional two 
points near the boundary are just inside the ellipse, more vor- 
tices will be present at large r± and these will lie at the other 
ends of the branch-cut lines; otherwise the branch-cut lines 
will extend to infinity. Fig. Ob) shows the case for 63 = 90°, 
giving a circle in parameter space. The number of vortices 
will therefore be close to the maximum for fixed values of r 2 
and r3 . 

The far-field solution is usually considered valid for a Fres- 
nel number Np = a 2 1 A$z ^ 1, where a is the largest trans- 
verse length scale present in the system (see, e.g., L29]). 
We take this to be the maximum pinhole-pinhole spac- 
ing. The far-field condition corresponding to Fig. is then 
(3Aq) 2 /Aqz <tc 1 or, say z > 100 Aq. The value of zo for 
numerical simulations in Fig. was deliberately chosen as 
25 Aq to show up visual disagreements between the vortices 
in the numerically determined phase map and their analyti- 
cally determined locations. Small arrows illustrate the true 
vortex positions for some vortices which do not coincide with 
their predicted far-field positions. The discrepancy between 
numerical and analytical results is highlighted in Fig. Oc)> a 
case arising from all parameter space lattice points lying close 



to the ellipse boundary. In Fig. Od), the Fresnel number is 
smaller and the correspondence is improved, but for any finite 
Zo, there will be some angle, close to n, at which the numer- 
ical and analytical prediction disagree by an arbitrarily large 
amount. Predictions are more reliable for (m, n) lattice points 
closer to the ellipse center. In contrast with the far-field pre- 
diction that vortices may abruptly appear and disappear with 
infinitesimal changes in source arrangement as lattice points 
cross the ellipse boundary, we observe through simulations at 
finite zo, that the vortex location becomes highly sensitive to 
source arrangement, with vortices arriving from and escaping 
transversely to an infinite r ± distance. This point is further 
explored in Sec. lIIIHl 

E. Source phase variation 

We now consider the effects of source phase variation on 
the geometric parameter space description. 

The functions M(m) and N(n) [Eqs (fTUb l contain the vortex 
coordinates' dependency on source phase via Eqs (IT3l) and 
(fTTb , where cpj is the relative phase of the 7th source and (pi = 
0. These may be rewritten as 

M(m) = 2tt [1 + 3 (m - (p 2 /2n)] (28a) 

and 

N(n) ee 2tt [2 + 3 (n - 3 /2tt)] , (28b) 

to highlight the 2n phase periodicity. 

We see that a 2n change may be absorbed as an integer 
change in an associated parameter- space coordinate m or n. 
The subtraction of a value from m or n corresponds geomet- 
rically to a translation of the (m, n) lattice points along the 
associated axis, with a change of 2n effecting a translation of 
one lattice unit. 

In summary, parameters defining the source configuration 
or the wave-number map to the parameter space as different 
ellipse constructions, whereas source phase variations corre- 
spond to translations of the lattice itself. 

F. Estimate of number of vortices 

In contrast to the case of three interfering plane waves re- 
viewed in Sec.HH in which infinitely many nodal lines are pro- 
duced, the analysis of the preceding sub-section implies that 
only a finite number of vortices are produced by three over- 
lapping spherical waves, provided that there is a finite spacing 
between the three corresponding point sources. Here, we give 
a simple means to estimate the number of nodal lines, as a 
function of the geometry of the three point sources. 

As mentioned earlier, the two-valued nature of Eq. (ITTI) im- 
plies that each (m, n) pair gives rise to a pair of vortices. Be- 
cause the m and n values have unit spacing [see Fig. Ob)], the 
number of vortices n v may be approximated by twice the area 
Tia'b' of the ellipse. Thus 

n v * 2ndb\ (29) 
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or 

k 2 I I k 2 £/ 
n v « — \r2r3 sin #3 = , (30) 

Z7T I I 7T 

where srf is the area of a triangle whose vertices coincide with 
the locations of the three point sources. 

One may ask whether it is possible to develop an exact 
expression for the number of lattice points enclosed by our 
parameter- space ellipse, thereby improving on the approxima- 
tion for n v given in Eq. (l30l) . Indeed, this question is addressed 
by a famous problem in number theory known as "Gauss's cir- 
cle problem". With 2 = 2n/3, 3 = -2n/3 and 3 = n/2, the 
conditions corresponding to the common form of Gauss's cir- 
cle problem are established; the ellipse is a circle and is cen- 
tered on a lattice point. The resulting solution for this case is 
both involved and well known, and so will not be given here 
(see, e.g., Andrews ll32ll ). Instead, we merely note that this 
establishes an unexpected and beautiful connection, between 
the physical system considered here, and a certain problem in 
the theory of numbers. 




FIG. 8: (Color online) Vortices move with varying source phase 
along trajectories and lie at the intersections of two independent sets 
of hyperbolas. Two independent phases give rise to two sets of hy- 
perbolas indexed by m (solid lines) and n (dashed lines) respectively. 
The hyperbolas are shown overlaying phases from Eq. (|6]). Note that 
the result of this equation is > for arbitrary (r ± , 6). The parameters 
are the same as Fig.^f). 



G. Vortex trajectories from phase variation 

In most of the preceding discussion, the phases <p\,<p2,<p3 
of the sources were all set to zero. Here, we consider how 
vortices move in response to varying the phase of one of the 
sources. 

The equations of these curves are found by eliminating the 
phase (pj corresponding to the source rj in Eqs (TT3T) and ([T71) . 
For example, the trajectories for variation of source r 2 are 
found by solving Eq. (TT3T) for 02 and substituting this into 
Eq. ([T71) . Repeating this for 3 gives a second set of trajec- 
tories along which the vortices move with 3 variation. The 
resulting equations are 



r ±m 



M(m) zo 



V(3£r 2 cos# m ) 2 -M(m) 2 



and 



r ±n 



N(n) zo 

V(3£r 3 cos(^-# 3 )) 2 -A^) 2 



(31a) 



(31b) 



where the subscripts mn have been changed to m or n to 
highlight the independence of the equations with respect to 
the complementary parameter space coordinate. Note that 
Eq. (13 lab is the same as Eq. ([T5K with minor manipulation. 

Fig. [8] shows the trajectories overlaying the numerical sim- 
ulation results for the source arrangement seen in Fig. Of). 



We have seen (Sec. IIII El) that variation of 02 (or 3 ) cor- 
responds to a translation of the m (or n) lattice alone. Thus, 
our choice to eliminate the respective phase has led to the sep- 
aration of the m and n dependence. The equations describe 
hyperbolas, where the polar coordinate system origin is cen- 
tred between the two hyperbola branches. Note that this polar 
form is less common than that typically seen in the study of 
conies, in which the origin is placed at the focus of one branch. 



Eqs (I3T1) may be thought of as a (one-to-many) mapping from 
lines of constant m or n in parameter space to hyperbolas in 
real space. 

The angle # 3 in Eq. J31bl) corresponds to the rotation of the 
^-indexed hyperbola axes with respect to the coordinate sys- 
tem axes. The square-root limits the range of real solutions, 
hence the number of hyperbolas and the number of vortices. 
The domain of the parameter- space variable m (or n) is re- 
stricted to match that determined by the parameter- space el- 
lipse. Thus these equations may be applied to determine the 
valid domains of m or n independently, without recourse to 
the ellipse solution. 

Fixing the source phases defines two particular hyperbo- 
las corresponding to the parameter space coordinates m and 
n respectively, with vortices located at the intersections of an 
m and an n curve. Both branches are included, giving four 
separate curves. Where one branch of an m (or n) hyperbola 
intersects both branches of an n (or m) hyperbola, only one 
intersection is observed to correspond to a physical vortex so- 
lution. 

These equations also describe an intersecting curvilinear 
coordinate system. The location of a vortex on the phase map 
may be thought of as addressable by selection of a particular 
patch having coordinates (m, n) and then finely addressable 
within that patch by variation of the source phases. Variation 
of the phases from to 2n shifts the vortices along one of the 
two families of curvilinear trajectories within the patch. When 
there are a moderate or large number of vortices present, the 
patches in the central region of the phase map are correspond- 
ingly small. In this case, the vortices shift a small distance 
with changing phase. In the outer, sparsely populated re- 
gions, small changes in phase lead to arbitrarily large distance 
changes since the outermost patches have boundaries at infin- 
ity. 

It is not unexpected that vortices should lie along intersect- 
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ing hyperbolas. Two point sources naturally give rise to sur- 
faces of constant phase that are hyperboloids of two sheets 
having the sources as their foci. The difference in path length 
being constant establishes the condition for constant phase 
and is equivalent to the well known geometric construction 
method for hyperboloids. The intersections of hyperboloids 
of two sheets with the z = Zo plane will be hyperbolas of two 
branches. These will not have phase values at the amplitude 
maxima or minima due to any individual pair, but will be at 
some other value x - a + ib on one set an d ~X - ~ a ~ ib on 
the other, summing to zero at the vortex locations. 



H. The parameter space near collinearity 

Here we examine the behaviour of the nodal lines as 63 ap- 
proaches 71 (three collinear sources) when 2 = 03 = 0, since 
we observe a marked deviation between numerical simula- 
tion and analytical results in this regime [see Fig.[6fc), which 
clearly illustrates this discrepancy] . 

An analysis around 63 = n may be performed by substitut- 
ing 63 = 71 + s into Eq. (fT2l) . where s is small: 



— ( cos(7r + s) + tan 6 sin(7r + s) | = M^—L m 
T2 \ ) Mini) 



(32) 



Applying the small-angle approximations cos s ~ 1 , sin s ~ s, 
Eq. OH) yields 



1 tr 2 N(n) 
trni6 mn (£) = - I — — : + 1 



s \ V3 M(m) 



(33) 



for the 6 mn coordinate of the vortex core. Substituting into 
Eq. dT31) and applying Eq. (IT6l) gives r ±mw for our special case 



zo 



1 



3kr 2 

M(m) ) 1 + tan 2 G m 



(34) 



- 1 



Finally, substituting Eq. (l33lh we get 



Zo 



(35) 



(3kr 2 /M(m)Y 



- 1 



s 2 W3 M(m) 



In the limit s — > 0, we now see that 6 mn (s) —> n/2 and 
r±mn(s) — > +izo- Regarding the latter limit, r ±mn (£) first ap- 
proaches infinity and then becomes imaginary. 

We may widen the context of this result by realizing that it 
is but one example of a lattice point crossing an ellipse bound- 
ary. Whenever this occurs, the far-field prediction is that as- 
sociated vortices will be created or destroyed instantaneously. 
In contrast, in numerical simulations, which are at some finite 
distance from the source, vortices rapidly enter from or escape 
to infinity. 



I. Relation to the Young's three-pinhole interferometer 

Here we show how the to map our results for three spheri- 
cal point sources, hitherto the main subject of this paper, onto 
a three-pinhole Young's interferometer. In this interferometer, 
coherent radiation illuminates a black screen that is punctured 
with three small pinholes, with the resulting transmitted radia- 
tion being observed at a distance that is large compared to the 
spacing between the pinholes. Note that the assumptions of 
equal amplitude and a single wavenumber, applied in Eq. ©, 
correspond to uniform illumination of the screen. The follow- 
ing argument is based on the Rayleigh-Sommerfeld diffrac- 
tion theory (see, e.g., B33I1 ). 

The Rayleigh-Sommerfeld diffraction integral of the first 
kind yields a rigorous solution to the time-independent 
Schrodinger equation (Helmholtz equation) in a vacuum-filled 
half-space, for a field that obeys the Sommerfeld radiation 
condition. For a field incident on an aperture A centered 
on the plane z = 0, the wavef unction U may be determined 
at an arbitrary point (x, y, z) in the half-space z > 0. With the 
boundary conditions that U(x,y,z = 0) « U^ l \x,y,z = 0) 
when (x,y,z = 0) is in A, and U(x,y,z = 0) « when 
(x, y, z = 0) is not in A, the diffraction integral reads: 



U(x. 
where 



;y, z ) = U (i) (x',y',0)K(x,y,x',y')dx'dy', (36) 



2n dz 



1 d ( exp (ikr) \ 



(37) 



is a propagator and r = -fix 
tion of the derivative gives 



■ x') 2 + (y - Y) 2 + z 2 . Evalua- 



ir lz ( ik ! 1 n ^ 
K= — exp (ikr) . 

Z7T r \ r r L ! 



(38) 



For a pinhole aperture, = S(x - x' ,y - y',z), and the 
wavefield U (x, y, z) takes the form of the propagator, Eq. (l38i 
In general, the pinhole aperture does not produce spherical 
waves. However, when the observation point is such that r » 
/io, the first term in parentheses dominates the second, giving 



K 



ikz 
2nr 2 



exp (ikr) . 



(39) 



It can now be seen that Eq. ©, which resulted from fac- 
toring out a common exp (ikr) jr term, would be unchanged 
by instead factoring out a common term of ikz exp (ikr) /2nr 2 . 
Similarly, when one is both in the far-field and close to the 
z-axis, z ~ r, so Eq. (|39t reduces to a spherical wavefunction 
with a constant multiplier, as asserted. 

In summary, the analysis of the preceding sections may be 
mapped onto the case of a Young's three-pinhole interferom- 
eter, in the far field. The nodal planes of the two-pinhole 
Young's experiment, which are unstable with respect to per- 
turbations, therefore decay into a nodal-line network of vortex 
cores when the third pinhole is added. 



10 



IV. CONCLUSION 

A network of phase vortices was seen to be generated by 
the superposition of three stationary- state sources of outgo- 
ing spherical waves. We presented an analysis of the structure 
of the associated vortex cores (nodal lines), in the far-field 
regime. A finite number of vortices was seen to be generated. 
Determination of the number of such vortices was mapped 
onto the problem of determining how many points, on a two- 
dimensional cubic lattice, lie within a given ellipse. The equa- 
tion of the ellipse depends in a known way on the geometry 
of the sources. The parameter space description also gives 
some insight into the effects of varying both the arrangement 
of the three sources, and their relative phase. Indeed, phase 
variation of two of the sources provides a means for precisely 
positioning one or several vortex cores. Lastly, we showed 



how to map all of the preceding analyses onto the problem of 
determining the far-field disturbance that results when a three- 
pinhole Young's interferometer is coherently illuminated. In 
contrast to the classical two-pinhole Young's interferometer, 
in which the resulting diffracted field vanishes over a series of 
nodal planes, the three-pinhole interferometer yields a quite 
different phase topology, permeated with a rich structure of 
nodal lines that thread vortex cores. 
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